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Abstract 

One may argue that QCD solves the strong CP problem by itself, without having to introduce 
new symmetries and particles. To test this idea, a lattice simulation is performed. The problem is 
investigated in the CP 3 model first. It is found that the model has a first order phase transition 
in 9 from a confining phase at small 9 to a deconfining phase at large 9, and that the critical value 
of 9 decreases towards zero as [3 is taken to infinity. This suggests that 9 is tuned to zero in the 
continuum limit. Preliminary studies of the SU(2) Yang-Mills theory in four dimensions show a 
phase transition in 9 as well, so that it is quite likely that the strong CP problem in QCD is solved 
along the same line. 



In QCD, as well as in other theories that possess 
instantons, the proper vacuum states are superpositions 
of vacua of different winding numbers n: 

|#) = ]Texp(k9n)|n},#e [0,2tt). (1) 

n 

These so-called 9 vacua are realized by adding a CP 
violating term to the action, 

S e = S - i9Q, (2) 

where S is the standard action and Q is the topological 
charge. A priori 9 is a free parameter. Since no CP 
violation has been observed in the strong interactions, 9 
must however be very close to zero, t This constitutes 
the strong CP problem. The present upper bound is 
9 < 10" 9 [1]. 

A popular picture of the QCD vacuum is that 
of a dual superconductor [2] in which color magnetic 
monopoles condense and color electric charges, i.e. 

f In the presence of quarks the 6 angle can be rotated into a phase 
of the quark mass matrix M and vice versa, due to the U^(l) 
anomaly. The effective angle is in this case 6 = 6 + argdetM. 



quarks and gluons, are confined by a dual Meissner 
effect. This picture has been successfully tested in 
lattice simulations [3]. In the 9 vacuum these monopoles 
acquire a color electric charge of the magnitude 9/2tt 
[4]. For 9^0 one would expect that the long-range 
color electric forces are screened by monopoles and that 
confinement is lost. So QCD would only be a viable 
(continuum) theory for 9 = 0. 

What makes a lattice simulation at non-vanishing 
values of 9 very difficult is the fact that the action is 
complex. As a result, standard lattice techniques are not 
immediately applicable. This has led us to investigate 
the problem in a simpler model first. 

A model, which in many respects is similar to QCD, 
is the CP^ -1 model in two space-time dimensions. 
The CP^ -1 model deals with N-component, complex 
scalar fields z a (x) of unit length: ~z a (x)z a (x) = 1, a = 
1, • • • , N. Out of these fields one constructs composite 
vector fields 



A p.( X ) = 7j z a(x) dp, z a (x). (3) 
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Figure 1. The free energy F(6) (a), the charge density q(0) (b) and the string tension <r(e, 6) (c) as a function of 
9 on the V = 64 2 lattice at (3 = 2.7. The solid curve in (a) is the prediction of the large-N expansion to leading 
order. Only the first half of the 6 interval is displayed. In the second half of the interval F(6) = F(2tt — 6). 



Then the action [5] can be written [6] 

S = f3 J d 2 x~D~, a {x)D ll z a {x), (4) 

where = + iA^. Thus the model describes a 
set of charged scalar fields interacting minimally with a 
composite gauge field. The topological charge is given 
by 

Q=7T ( d 2 xc^d^A„(x) = -!- / d 2 xF 01 . (5) 
Zir J Zir J 

The 9 dependence of the theory is governed by the 
partition function 

Z(ff) = Y, e M^Q)p(Q) = eM-VF(6)), (6) 
Q 

where p(Q) is the probability of finding a field 
configuration with charge Q, and F(9) is the free energy 
per space-time volume V . In terms of F(9) the average 
topological charge density is given by 

l(e\Q\e) = -i q (e) = -i^l, (7) 

and the string tension of two external particles of charge 
e and — e (in units of the intrinsic charge) turns out to 
be 

a(e,9) = F(9 + 2ire) - F{6). (8) 

We have chosen to investigate the CP 3 model. For 
details of the calculation see ref. [7, 8]. 



If our idea is correct, we should find a first order 
phase transition in 9 from a confining phase to a Higgs 
or Coulomb phase. On a finite lattice and at a finite 
value of [3 the phase transition is expected to occur at 
a value 9 = 9 c (p,V) > 0, where V now is the lattice 
volume. Only on an infinite lattice and at [3 = oo, i.e. 
in the continuum limit, would we expect that 9 C = 0. 

A first order phase transition manifests itself in a 
kink in the free energy, as well as in a discontinuity in 
the first derivative of the free energy, i.e. the topological 
charge density. In Fig. 1 I show F(9) on the V = 64 2 
lattice at [3 = 2.7. For comparison I also show the 
prediction of the large-N expansion to leading order 
[6], F(9) = const. • 9 2 . We see a distinctly marked 
kink at 9 = 9 C « 0.5 tt : while F(9) increases roughly 
proportional to 9 2 up to 9 = 9 C , F(9) turns out to be 
constant (within the error bars) for 9 > 9 C . In Fig. 1 
I also show the topological charge density on the same 
lattice. According to (6) this can be interpreted as a 
background electric field. We see that q(9) increases 
almost linearly with 9 up to 9 = 9 C , where it jumps to 
zero and then stays zero over the rest of the interval 
(again within the error bars). Thus the phase transition 
is marked by a collapse of the background electric field, 
presumably due to pair production. 

To show that the phase transition is a deconfining 
phase transition, I have plotted the string tension in Fig. 
1, again on the 64 2 lattice at (3 = 2.7, for two different 
charges. To keep the discussion simple, I have restricted 
myself to the case of small fractional charges. This 
matters because the vacuum will change its properties 
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Figure 2. The critical value C (P, V) as a function of V 1 for 
three values of f3. The lines are a linear fit to the data points. 

if it is exposed to a strong external electric field. We see 
that the string tension is zero (within the error bars) for 
6 > 6 C . This result is a consequence of the property that 
F(6) is constant for 6 > 6 C . The string tension reaches 
its peak at 6 = 6 C — 2ire. 

In Fig. 2 I show all our results for 9 C {P, V) for 
a variety of lattice volumes ranging from V = 28 2 to 
V = 200 2 and for three values of /?: /? = 2.5, /? = 2.7 
and /? = 2.9. The lattice data display strong finite size 
effects. For a first order phase transition we expect 

6 c (p,V)-6 c (p,oo)<xV- 1 (9) 

for fixed (3. When 6 c (j3,V) is plotted against V -1 , we 
find that for all three values of [3 our data fall on a 
straight line, in accordance with the predictions of a first 
order phase transition. This allows us to extrapolate the 
lattice results to the infinite volume. 

We find that the extrapolated values of 6 C fall rapidly 
to zero as [3 is taken to infinity. Our results are 
consistent with a decay like 6 C oc where £ is the 

correlation length. (For our values of [3 the correlation 
length is 4.5, 8.8 and 18.5, respectively.) In the strong 
coupling limit the model is known to have a first order 
phase transition at 9 = 6 C = 7r [9] . This suggests a phase 
diagram of the form shown in Fig. 3. 

Thus 6 = is the only value at which one can take 
the continuum limit, at least within the confining phase. 
This resolves the strong CP problem. 

As a next step towards solving QCD, we are 
currently investigating the problem in the SU(2) Yang- 
Mills theory in four dimensions [10]. Preliminary results 
indicate a first order phase transition in 9, just as in 
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Figure 3. The phase diagram. The horizontal line at 6 = 7r is 
the strong coupling prediction, and the symbols are the lattice 
data extrapolated to infinite volume. 

the case of the CP 3 model. Because of lack of space 
I cannot show a figure here. To establish the phase 
diagram firmly, we will have to repeat the calculation 
on lattices of various sizes and for several values of [3 
still, as we did for the CP 3 model. 
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